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Abstract 

Inspired by scientific collaboration networks, especially our empir- 
ical analysis of the network of econophysicists, an evolutionary model 
for weighted networks is proposed. Both degree-driven and weight- 
driven models are considered. Compared with the BA model and other 
evolving models with preferential attachment, there are two significant 
generalizations. First, besides the new vertex added in at every time 
step, old vertices can also attempt to build up new links, or to recon- 
nect the existing links. The reconnection between both new-old and 
old-old nodes are recorded and the connecting times on every link is 
converted into the weight of the link. This provides a natural way for 
the evolution of edge weight. Second, besides degree and the weight of 
vertices, a path-related local information is also used as a reference in 
the preferential attachment. The path-related preferential attachment 
mechanism significantly increases the clustering coefficient of the net- 
work. The model shows the scale-free phenomena in degree and weight 
distribution. It also gives well qualitatively consistent behavior with 
the empirical results. 
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1 Introduction 



Network analysis is now widely used in many fields ^ [2] • Recently more 
and more works on weighted networks appears in both empirical and mod- 
elling analysis. In a weighted network, the weight on the edges provides a 
natural way to take into account the interaction strength, while in a binary 
network, the edges only represents the presence or absence of interaction. 
This capability will probably carry more information about the interaction. 

The first problem before any analysis can be applied to the weighted 
networks is how to assign the weight to edges. This problem is quite non- 
trivial. Several ways to assign the weight have been introduced. One is 
to transfer some quantities from non-weighted networks into the weight of 
edges. They are usually related to the degree or other intrinsic quantities of 
the nodes 0111 Ej- As in 0, the weight of an edge is measured by the point 
degree of its two ends, which are defined phenomenologically from binary 
networks. It is helpful to describe new properties of the binary networks, 
but it includes no more information than the origin binary networks. Some- 
times, the real- world phenomena investigated provide a typically natural 
measurement of the weight, such as the number of flights or seats between 
any two cities in airport networks [HI 13 IB1> the reaction rate in metabolic 
networks [H] and so on. In the works of modelling weighted networks, weights 
on edges are generated from priori distribution jlUl ITT1 112j . From the view 
point of empirical study, we never know such models already acquire the 
real structure of weighted networks or not. 

However, some weighted networks such as scientific collaboration net- 
works are different with the above networks. In the collaboration networks, 
the connection times is a natural quantity which is related very closely to 
the weight. But there is no explicitly expression between this quantity and 
weight. Let's think about the times of coauthoring between two scientists. 
Obviously, more times represents closer relationship in the sense of trans- 
portation of scientific ideas. Therefore, in scientific collaboration networks, 
usually the happening times of the event is converted as the weight of the 
edge. Yet different authors may use different expressions [7| IT5| I14| 1151 H5j . 
As to which definition behaviors better, and whether or not there are some 
general rules to define weight, we do not have the final answer yet. 

The second problem related to the weighted networks is how to extract 
information from weighted networks constructed by the above ways. Es- 
pecially one may concern about what's the role of weight, or what's the 
significant difference brought by weights compared with binary networks. 
In order to answer the above questions raised from those two aspects, we 
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have to consider the third problem, modelling the weighted networks. 

For instance, by investigation of modelling works, if we find that in 
order to construct a well-behavior model of weighted networks, the degree is 
the only variable directly coupled with evolution, while the weight is never 
needed to directly be brought into the evolutionary process, then we may 
think that the weight just relies on a higher level structure. The weight is not 
crucial in network analysis although it is important and necessary. Or quite 
the contrary, if in order to get weighted network behaviors consistently with 
real phenomena in the modelling work, the weight must be coupled directly 
with the evolution. Then the weight should play a significant role in the 
way to extract information from weighted networks. 

In this paper, we tried both degree-referred preferential attachment and 
weight-referred preferential attachment in our evolving model and compared 
the results with the empirical analysis from |141 1151 116j . 

Actually, there are already many evolving models for weighted networks. 
Some models introduced prior weights into edges with the evolution of net- 
works. In each link j «-> i from the newly added node j is assigned a 
weight as uiji = fc ., ; where {i'} represents a sum over the m existing 

nodes to which the new node j is connected. Zheng[l] has improved this 
idea. In his model, the weight of a link depends not only on total degree 
of the existing nodes, but also on some intrinsic quality ("fittness") of the 
nodes. In ^7j, the weight of a link depends on randomly modified intervals 
between the time at which linked vertices are connected to the system. In 
[P2] . the weight Wij of a link Zy connecting a pair of nodes (i and j) is defined 
as Wij = (wi + Wj)/2, and W{ is defined as i node's assigned number (from 
1 to N) divided by N. In some evolving models [TU1 lllj . the weight w is 
assigned to the link when it is created and it is drawn from a certain distri- 
bution. As pointed out in [T3], most models here are not really evolutionary 
models in the sense of weight. The weight keeps the same value after it was 
assigned onto its edge. Or some extra quantities are introduced to drive the 
evolution of networks. 

Recently, some evolving models are set up in which the weights are cou- 
pled directly with the network evolution. In the paper ^j, a weight-driven 
model was proposed and the weight of link changes with the network evo- 
lution. In this paper, the new edges starting from the new vertex added 
in at every time step are preferentially attached to old vertices determined 
by their strength, or vertex weight. After the attachment, an increase of 
weight 5 is distributed among all the edges connecting to the chosen old ver- 
tices. The model yields a nontrivial time evolution of vertices' properties and 
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scale- free behavior for the weight, strength, and degree distributions. In the 
paper 19 , Bianconi has presented a model with co-evolution of link weight 
and strength. In his weighted fitness network model, the fitness of node 
and link and results in the structural phase transition of the network are 
introduced. In [2U|, the network evolves with connectivity-driven topology 
and with the weight assigned from a special distribution pk{x) of weights. 

Although the models mentioned above coupled the weight and network 
evolution, we think that the dynamical process of the weight in Barrat's 
model [TH] is quite artificial or say not very general, or like in the other two, 
extra quantities not rooted in network has to be used. The authors of |18j 
gave some arguments for this as to justify the process from the background 
of airport network [TH]. But they took weight as a quantity independent 
on connecting. However, as we have mentioned before, weight usually re- 
lated closely to connecting times. Especially for the actors and scientists 
collaboration networks, using weight converted from connecting times is a 
convenient way to construct weighted networks. Therefore, it seems that 
such a pure weight-driven model depends too much on this artificial dy- 
namical process of weight. Now, our empirical investigation on scientific 
collaboration networks give us some hints on modelling weighted networks. 

In our model, we keep the relationship between weight and connecting 
times, and only quantities directly rooted in networks are used. So the 
picture of the evolution looks like the connecting times evolve according to 
weight, and then the new connecting times comes into the weight, which 
drives the evolution of the system again. Or in our degree-driven model, 
connecting times evolve according to degree, and degree increase due to 
connecting, and then all the connecting times are recorded and converted 
into weight. 

Another important improvement of our model is the introduction of 
local-path-related preferential attachment, the 5 term in our model. This 
mechanism works for the network evolution in the real world but is ne- 
glected by other models. It is helpful to increase the clustering coefficient 
of the networks. One major difference between empirical results and most 
models is about the clustering coefficient. Usually, BA model 2 1 or simi- 
lar models ^H]) given a quite low clustering coefficient while in reality, real 
phenomena show highly clustered behavior. Of course, the WS model of 
small world network |22j gives high clustering coefficient because it starts 
from a regular network, not on the way of evolutionary network models. 
Some evolutionary models do give high clustering coefficient |23l l2*H I27j . In 
[2*3] , if an edge between v and w was added, then add one more edge from 
v to a randomly chosen neighbor of w. In [21], one randomly chosen per- 
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son introduces two random acquaintances to each another who haven't met 
before. Another idea is to introduce an extra Euclidean distance, and ver- 
tices prefer to interact with nearby vertices. Therefore, in order to increase 
the clustering coefficient, new mechanisms which are not rooted directly in 
the network have been introduced. Now, we introduce the 8 mechanism, 
where all quantities still comes directly from the topological structure. This 
requires no more extra information, but just a little knowledge about the 
local structure. Here 'a little' means one only need to know the information 
about the second, or third nearest neighbors, not any more. 

The detailed comparison will be done between the results from the mod- 
els and our empirical results from |15| I16j . The description of the general 
model is given in Section SJH The asymptotic distributions of vertex weights 
for the weight-driven case is also given analytically in Section and they 
are well consistent with results of numerical simulations. In Section [JIJ in 
order to compare with the empirical study of econophysists collaboration 
network, we extend our model onto directed weighted networks. In this 
comparison, they show nice agreement. 

2 Models and theoretical analysis 
2.1 The model 

A N- vertex weighted network is defined by a N x N matrix ujy, which 
represents the weight on the edge from vertex i to j. Similarity weight is 
used here. So the larger the weight is, the closer the relation between the 
two ends nodes are. wi m = means no relation between vertex I and m. 
Suppose the edge weight Wij is related to the connecting times 2y between 
vertex i and j, by 



such as the tanh function iUy = tanh (aTy) we used in ^3E3> or just linear 
relation = aTij used by other authors [71 HI]. 

Our most general model is given as follow. Starting from a fully con- 
nected no initial network, with initial times Tjj = 1 (and initial weight 
Wij = f (1)), at every time step, 

1. One new vertex is added into this network, and I old vertices are 
randomly chosen from the existing network. 

2. Every one (denoted as vertex n) of them can initially activate a temp- 
tation to build up m connections. The probability for every link from 





5 



n connecting onto vertex i is given by 



IU^i = (l-p)^ r + ip-S)=^- + 6 lm , (2) 

where ki is the degree of vertex i, Wi = Y2j w ji is the "onto" vertex 
weight of vertex i, Z n j is the similarity distance from n to i, and 
d% means the dth neighbors of vertex n. For example, d\ is the set of 
nearest neighbors, d\ means the second nearest neighbors, so that dn 2 
in the expression refers to both of them. Intuitively, similarity distance 
means the maximum distance between two vertices because the weight 
is defined here in the way that the larger the closer. Usually, in calcu- 
lation of network analysis, the dissimilarity distance corresponding to 
the shortest distance is used more often. 

3. After we got an end node i* chosen from all vertices over the exist- 
ing network by the probability above in equ(J2), the connecting times 
between vertex n and i* increased by 

T ni *(t + l)=T ni *(t) + l. (3) 

4. The weight of the edges changes as 

w nt *(t + l) = f(T ni *(t+l)). (4) 

Although our general model defined above can be applied to directed net- 
works, in the following analysis we assume that Wij = Wji. An increase on 
Tij immediately reflects another increase on Tji. Except for the comparison 
with empirical results, on most cases, the linear function is used for the 
relationship between connecting times and weight for the simplicity, 

Wij = aTij. (5) 

2.2 Analytic results of the weight distribution 

Now we try to get the analytical results for the vertex weight distribution 
under the simplest weight-driven model. For the link weight given by equ([5|) . 
the weight of vertex is given by 

3 
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we suppose that the newly added vertex and the old vertices are informed 
of the weight of the other vertices and the network is pure weight-driven. 
In this case it is attached with preferential linking described by p = 1 and 
5 = in equ©, that is the the connection probability is 

n_ = =B- (7) 

The master equation for the evolution of the average number of vertices with 
weight w at time t is 

N(w,t + 1) =N( W ,t) + m • (1 + • {w - iyN ^ W l%^ {w ' t] . 

l-N(w,t) l-N(w-m,t) W , c 
N ' N ' " w > m 

Here Y,w w ■ N ( w > t - !) = 2 • E o + 2 ■ m ■ (! + ' 1 is the total weight 
and = uq + t is the size of system at time t. The equation describes 
the increasing of preferential linking since the new vertex is added and old 
vertices are selected. The first term reflects the preferential attachment 
(|T|) used to select the other end of the link, while the following two items 
correspond to the random selection of I old vertices. When Eq and no are 
much smaller than t, the size of the network N is approximately the time 
steps t. Then the master equation® can be written as 

(i + 1) • p(w, t + l)=t- p(w, t) + ±-[(w-l)- p(w - l,t) 

(9) 

-k ■ p(w, t)]-l ■ p{w, t) + 1 ■ p(w -m,t) + 6 Wjm 

where p(w, t) ~ N (™& \ s the density of vertices with strength w at time 
When t is larger(i S> 1) enough, 

(t + 1) • p(w, t + l)-t- p(w, t) = p(w, t). (10) 

We get from equ© 

— — — — — I ■ [p(w) — p(w — 1) + p(w — 1) — pi ic — 2) 



2dw 

+p{w — 2) — • • • + p(w — m + 1) — p(w — m)] + 5 U 



-d(wp(w)) , dp(w) , r 



(11) 
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Figure 1: Comparison between the numerical solution results from eq© 
and the analytical results of ea ()14j) for different I. Other parameters are 
tiq = 10, m = 5. 



This is further written as 

d A 



p(w) + —— [- • w ■ p(w) + I ■ m ■ p(w)] = 5 V 
aw 2 



For w / m, we get 



(2-l-m + w)- = -3-p(w) 

aw 



We arrive at the final vertex weight distribution 

/)(((•) <x (w + 2 ■ I ■ • 



m) 



(12) 



(13) 



(14) 



In Fig(fJQ), we compare the numerical solution of equ© with the analytical 
results eau (|14|) . it shows a nice consistence. We can find that the lower end 
of strength is obviously affected by the parameter I and m and departure 
from power-law, while the upper end is still distributed as power-law. In the 
section ^3.21 we will compare the analytical results with that of computer 
simulation in Fig(??). They are also consistent very well. 



3 Numerical results 
3.1 Degree-driven Model 

First, we consider a special case, p = and 6 = 0. In this case, our model 
is fairy similar with BA model, except now, besides the new vertex added 
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in, the old vertices can also be activated. This assumption has been used 
already in several evolving models especially for the modelling of cooperation 
networks |26l 1271 128j . Another difference between this case and BA model is 
that the reconnection of link is allowed and recorded. Later on, it will be 
converted as the weight of link. Fig(J2J shows the typical behavior of degree 
and vertex weight distribution. The weight distribution of links obeys also 
power law as shown in Fig©. These results are consistent with the typical 
result from empirical studies qualitatively, such as distribution of vertex 
weight for airline networks [HI El and collaboration networks 7 . 

Compared with BA model, the introduction of parameter I is new, so 
how I will effect the behavior of the model? One limit situation is when 
/ = 0. All the contributions to the weight come only from the new ver- 
tices. So our model comes back to BA model except some new links may 
be repeated. The degree distribution is the same as BA model. The vertex 
weight distribution is almost the same as degree distribution but there is no 
power-law distribution of edge weight at all. The increasing of I will affect 
the degree distribution. The lower end will departure from the power-law 
distribution but show the "droop head" shape observed in many empirical 
studies. Another limit situation is / ^> 1. In that case, the increase of in- 
ternal links has much more effects on the network evolution compared with 
the growth of the network. The network will lose the power-law behavior in 
the lower end, although in a quite large domain of I, the power-law behavior 
of degree and weight distributions are robust, especially in the upper end. 
The effect of I is shown in FigQ. 

All the results are the average of 10 simulations for different realization 
of networks under the same set of parameters. The network sizes are all 
reach 10000 nodes. We have compared the distribution with that of the 
network with 50000 nodes. They are almost the same so that a network 
with 10000 nodes can give us a nice description for asymptotic distribution. 

3.2 Purely weight-driven model 

Now we assume that the vertex weight plays the most significant role in 
the evolution so that p = 1 and 5 = 0. This means the scientists choose 
their cooperators according to the weight, instead of focusing on degree. 
Therefore, weight is the fundamental character of vertices. Intuitively, the 
meaning of the degree looks like the extensiveness of the working style while 
the weight considering both extensiveness and intensiveness. So it's not very 
surprised that weight can unconsciously be used as a scale to attract more 
cooperators. In fact, the same idea of this weight-driven mechanism has 
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Figure 2: Degree (a) and vertex weight (b) distribution from degree-driven 
model. They are both power-law distribution with almost the same expo- 
nent, 7 = —2.7. The difference between these two distributions in upper tail 
is shown in (c) by Zipf plot. But most points are in the lower region. In 
this simulation, no = 10, m = 5, / = 1. 
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Figure 3: Edge weight distribution. Empirical results from Newman's data 
of scientific collaboration networks (a) and collaboration network of Econo- 
physists (b). Simulation result of the model is shown in (c). 
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Figure 4: (a) Degree distribution affected by the parameter I. (b)Larger 
I (I = 10, 40) makes the power-law behavior lose at the lower end. Other 
parameters are no = 10, m = 5. 

been used in Barrat's paper^HJ. As we mentioned in introduction, the only 
difference between Barrat's model and this model is the evolution mechanism 
of weight. In Barrat's model, it evolves directly by a phenomenological 
rule as a 5-increase, while in our model, it evolves indirectly through the 
connecting times T^j. 

The numerical results are given in Fig(J5J). They are consistent with the 
theoretical analysis. It has been found that those two limit cases seems 
have the similar qualitative behavior. However, as we mentioned in the 
introduction, the difference between those two cases, and the conclusion 
that which one behaviors better, implies the answer to the question that 
which quantity is the more fundamental one between degree and weight. 
Or put it in another way, should weight be a high-level quantity defined 
by degree, betweenness, whatever the basic network quantities, or directly 
from event represented by the network? Therefore, an conclusion about this 
comparison is essential for this issue. However, so far those two models under 
the limit cases provide the similar behavior. In the next section £QJ when 
the models are extended onto multilevel relationships to do a comparison 
between models and empirical results, at first we extended both those two 
limit cases. Both of them provide consistent behavior with the empirical 
results. After the similar results are found, only weight-driven model is 
compared with empirical results further. Of course, one can even try to 
study the behavior for a general p besides those two limit cases. However, 
since we have no way to make significant difference between the special 
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Figure 5: Vertex weight distribution for the purely weight-driven model: 
(a) comparison of the simulation results with analytical results; (b) the effects 
of different values of I. Other parameters are no = 10, m = 5. 




Figure 6: Comparison between purely degree-driven and purely weight- 
driven models on weight distribution. We could see that those two models 
behavior very similarly (a) except for the upper tail in Zipf plot (b). The 
parameters for two simulations are the same: no = 10, m = 5, I = 1. 
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cases, no further investigation about the general model has been done in 
this paper. To compare those two models with more empirical results and 
find the significant difference between them is really valuable for network 
analysis. 

3.3 Clustering structure when 5^0 

The mechanism represented by the 5 term in our model has taken the local 
information of the network into account. When a connection is built up by 
the active vertex, the probability of a vertex being chosen as the end vertex 
is higher if it has the closer relation with the attempting one. Therefore, 
hopefully, this will increase the clustering coefficient. In the original BA 
model, this mechanism is neglected. So it's interesting to just keep the 
degree term and the 5 in our model, and explore whether such mechanics 
increase clustering coefficient or not. 

3.3.1 5-mechanism applied onto model of non- weight networks 

For the non- weighted networks, every edge has the weight 1. The l n i in 
5 term is the similarity distance, which is the reciprocal of the shortest 
distance between vertex n and i. For the first link from new attempting 
vertices, just the BA rule of preferential attachment is applied, but for the 
links afterwards and the links from old vertex, the end point is determined 
preferential by both its degree and closeness with starting point. In Fig(jJJ), 
the simulations show this 5 mechanism significantly increases the clustering 
coefficient while the power-law distribution of degree still holds. 

3.3.2 5- mechanism in weighted network 

For the weighted networks, l n i in 5 term is the similarity distance as men- 
tioned before. In this case, just for simplicity, we consider the purely weight- 
driven model, which means p = 1, 5 ^ in equ©. Its effects on clustering 
coefficient are also shown in Fig0. The clustering coefficient under this 
mechanism reaches a stable value after a period of evolution. We can see 
that under some value of 5, for example, S = 0.8, the clustering coefficient is 
near 0.25, which is much higher than BA model, and even comparable with 
empirical results [T^l ITol |2*7| and results from other models 231 l2*ij. 

Since we have already shown that the 5-mechanism can increase the 
clustering coefficient, from now on, we will focus only on the comparison 
between simulation and empirical results on the distribution of degree and 
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Figure 7: Clustering coefficient for degree-driven (a) and weight-driven (b) 
model. Parameters are no = 10, m = 5, I = 1. 

weight, instead of on the clustering coefficient. So in the following discussion, 
5 is set to be again. 

4 Extended model and comparison with empirical 
results 

In the real world, relations of nodes usually are more than one levels and 
different relations have different contributions to the weight of link. For in- 
stance, in the empirical analysis in |15l I16| . we consider both co-authorship 
and citation as the ways of scientific idea transportation with different con- 
tributions. Even the worse is that citation is a directed network. So in order 
to compare the results from our models with the empirical studies, we must 
extend our current model into a multilevel directed networks model. There 
are two kinds of connecting times T^, where u = 1, 2 refers to co-authorship 
and citation respectively. Here the relation between connecting times and 
the link weight is given by a tanh function. The reason we prefer the tanh 
function in empirical studies is that, first, it has the saturation effect, which 
makes the contribution less for larger connecting times; second, it normalizes 
the maximum value to 1, which is the usual strength of edge in non- weight 
networks. So the two T^ - are converted into a single weight by 

Wi i = \^Z tanh ( a ^i) ' ( 15 ) 
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so that Wij is normalized to 1. And the probability distribution to chose the 
end vertex is consistently transformed as 




(16) 



while YlaP 11 = 1- Or in degree-driven model, 



n 



rt 



■i — 



5> 



(17) 



After vertex i* are chosen as the end vertex of a relation u between n and 
i* according to above probability distribution, the connecting time evolutes 
as 



For = 1, after that we need to set Tf* n (i + 1) = T^* (t + 1). For \i = 2, we 
skip this step. From its definition eau(|15|). the weight here is an integrated 
variable. This implies those two events can be triggered by each other, not 
developed separately. 

As it will be shown in Fig©, we have not found any significant differ- 
ence between those two models, so later on, when we compare models with 
empirical results, only weight-driven model are used there. As explained in 
the introduction section ^ measuring the role of weight by evolutionary 
models is one of the goals of our research which has not been achieved so 
far. We hope more comparisons between the behaviors of those two models 
and more empirical results will give an conclusive answer for this question. 

For directed network, the degree is divided into three quantities: out 
degree, in degree and total degree. For example, the in degree kf 1 is the 
sum of edges ending at vertex i, that is kf 1 = Y^j sign (wji). The out degree 
and weight are calculated similarly and the total degree and weight are the 
sum of in and out. The same situation happens to vertex weight, so there 
are out weight, in weight and total weight of vertices. From the simulation 
results, we can find that the total, in, and out degree and weight are all of 
power law distribution, as shown in Fig©. 

The more important comparison we want to do is between simulation 
and empirical results, especially on link and vertex betweenness, because 
they are global properties related with the whole structure of the networks. 
The link and vertex betweenness and their distribution could be gotten from 
the set of effective pathes between any two nodes. The average distance d 




(18) 
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(c) k (d) w 

Figure 8: Typical distributions of degree and weight (the three curves are 
In, Out and Total respectively) from pure degree-driven model (a), (b) and 
pure weight-driven model (c), (d). Again, the pure weight-driven model 
gives very similar results. 
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is defined as before as, 



ij 

in which, dij is the similarity distance of an effective path between vertex 
the larger the closer and equals to if no path exists. In fact, the 
above formula is not exactly the same as the one for non-weighted networks. 
First, because of the direction of edges, the number of total edges are now 
N (N — 1) instead of N ^~ 1 > ■ Second, the algorithm to search for such dij 
is slightly different with the usual shortest path in non- weighted networks. 
One way to make use of the shortest path algorithm is to transform the 
similarity into dissimilarity weight, so that the shorter the closer, and then 
use the usual shortest path algorithm to find all the distance. After that, 
transform it back into similarity distance. However, this is just an algorithm 
problem, has nothing to do with the structure analysis. 

In order to check the model, we compare the results with empirical re- 
sults from Econophysicists network, which has mostly been given in |151 116j . 
We compare the distribution of quantities of Econophysicists networks with 
numerical simulations, such as degree, vertex weight, vertex betweenness 
and link betweenness. It is interesting that the results are consistent well. 
It seems that the model reveals some basic mechanisms of the evolution of 
collaboration network. The parameters we used here for this comparison 
are I = l,p l = 0.2, p 2 = 0.8. 

5 Concluding Remarks 

In this paper, we presented an evolutionary model for weighted network, 
which integrates the contributions from both new vertices and old vertices. 
The two mechanisms, degree-driven and weight-driven preferential attach- 
ment are discussed, and both show a good consistence with empirical re- 
sults from network of econophysicists. Also a new mechanism, named as 
local-path-related preferential attachment, which makes use of some locality 
information is introduced here to increase the clustering coefficient of the 
network. Weight has been assigned to each link according the connecting 
times of the link, so the weight of link changes as network evolutes. Includ- 
ing the behaviors from the old vertices, recording all the connecting times 
and converting them into weight, are the most essential steps in our model. 
The way to incorporate locality information into network evolution by the 
(5-mechanism is also one point of this paper. 
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Figure 9: Zipf plot of degree and weight from empirical studies and sim- 
ulations, (a) and (c) are the empirical results of degree and weight distri- 
bution, (b) and (d) are the simulation results of the model. The model is 
simulated under the parameters: no = 10, m = 5,/ = l,p = 0.2, p 2 = 0.8. 
a 1 = 0.7, a 2 = 0.3 in eqi il5l The size of simulated network is N = 819 that 
is the same as empirical studies. 
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(a) Rank of Link (b) Rank of Vertex 

Figure 10: Zipf plot of betweenness for empirical studies and simulations, 
(a) Link betweenness. (b) Vertex betweenness. 

However, one of the most important tasks of this paper, which is to 
determine the role of weight and comparing it with the role of degree, has 
not been done yet. Although the comparison so far could not distinguish the 
degree-driven model and weight-driven model, we hope further comparison 
with empirical results will give a conclusive answer for this question. 
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